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Abstract 



A search for microscopic black hole production and decay in pp collisions at a center- 
of-mass energy of 7 TeV has been conducted by the CMS Collaboration at the LHC, 
using a data sample corresponding to an integrated luminosity of 35 pb^^. Events 
with large total transverse energy are analyzed for the presence of multiple high- 
energy jets, leptons, and photons, typical of a signal expected from a microscopic 
black hole. Good agreement with the expected standard model backgrounds, domi- 
nated by QCD multijet production, is observed for various final-state multiplicities. 
Limits on the minimum black hole mass are set, in the range 3.5 - 4.5 TeV, for a variety 
of parameters in a model with large extra dimensions, along with model-independent 
limits on new physics in these final states. These are the first direct limits on black hole 
production at a particle accelerator. 
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One of the exciting predictions of theoretical models with extra spatial dimensions and low- 
scale quantum gravity is the possibility of copious production of microscopic black holes in 
particle collisions at the CERN Large Hadron Collider (LHC) [1. 2J. Models with low-scale 
gravity are aimed at solving the hierarchy problem, the puzzlingly large difference between 
the electroweak and Planck scales. 

In this Letter we focus on microscopic black hole production in a model with large, flat, ex- 
tra spatial dimensions, proposed by Arkani-Hamed, Dimopoulos, and Dvali, and referred to 
as the ADD model |3, 4J. This model alleviates the hierarchy problem by introducing n extra 
dimensions in space, compactified on an n-dimensional torus or sphere with radius r. The mul- 
tidimensional space-time is only open to the gravitational interaction, while the gauge interac- 
tions are localized on the 3 + 1 space-time membrane. As a result, the gravitational coupling 
is enhanced at distances smaller than r, and Newton's law of gravitation is modified at short 
distances. The "true" Planck scale in 4 + n dimensions (Md) is consequently lowered to the 
electroweak scale, much smaller than the apparent Planck scale of Mpi ~ 10^^ TeV seen by a 
3 + 1 space-time observer. The relationship between Mo and Mpi follows from Gauss's law 
and is given as Mp^ = S/rM^+^r", using the Particle Data Group (PDG) definition |5|. 

Such a change in space-time structure and subsequent strengthening of the gravitational field in 
the ADD model could allow black hole formation in particle collisions at energies greater than 
Md, rather than Mpi, which is the case for a truly 4-dimensional world. Colliding particles 
would collapse in a black hole if their impact parameter were smaller than the Schwarzschild 
radius of a black hole with the mass Mbh equal to the total energy accessible in the collision. 
The Schwarzschild radius of a black hole with mass Mbh embedded in 4 + n space-time can be 
found by solving Einstein's general relativity equations and is given by miZl: 
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The parton-level cross section of black hole production is derived from geometrical considera- 
tions and is given by c ~ zrrg [1 , 2J. At LHC energies, this cross section can reach 100 pb for Md 
of 1 TeV. The exact cross section cannot be calculated without knowledge of the underlying the- 
ory of quantum gravity and is subject to significant uncertainty. It is commonly accepted HUll 
that the minimum black hole mass cannot be smaller than Md; although the formation 
threshold can be significantly larger than this. When a black hole is formed, some fraction of 
the colliding parton energy may not be trapped within the event horizon and will be emitted in 
the form of gravitational shock waves, which results in energy, momentum, and angular mo- 
mentum loss. This effect is particularly model-dependent for black hole masses close to Md. 
In general, black holes in particle collisions are produced with non-zero angular momentum, 
which also affects their properties and production cross section. 

Once produced, the microscopic black holes would decay thermally via Hawking radiation 
im, democratically (with equal probabilities) to all standard model (SM) degrees of freedom. 
Quarks and gluons are the dominant particles produced in the black hole evaporation (~ 75%) 
because they have a large number of color degrees of freedom. The remaining fraction is ac- 
counted for by leptons, W and Z bosons, photons, and possibly Higgs bosons. Emission of 
gravitons by a black hole in the bulk space is generally expected to be suppressed [H. In 
some models the evaporation is terminated earlier, when the black hole mass reaches Md, 
with the formation of a stable non-interacting and non-accreting remnant. The Hawking tem- 
perature for a black hole in 4 + n space-time is given by HHHIZl: Th = (in Planck units 
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h = c = ks = 1, where kg is the Boltzmann constant) and is typically in the range of a few 
hundred GeV. The lifetime for such a microscopic black hole is ~ 10^'^^ s |[lJ|2H7l. 

Here we consider semi-classical black holes, whose properties are similar to those for classical 
black holes described by general relativity and whose mass is close enough to Md so that 
quantum effects can not be ignored completely. There are also models ||T0l - [T2l of quantum black 
holes that decay before they thermalize, mainly into two-jet final states. We do not consider this 
signature here, leaving it for dedicated searches in the dijet channel ||T3l[T4| . 

The microscopic black holes produced at the LHC would be distinguished by high multiplicity, 
democratic, and highly isotropic decays with the final-state particles carrying hundreds of GeV 
of energy. Most of these particles would be reconstructed as jets of hadrons. Observation of 
such spectacular signatures would provide direct information on the nature of black holes as 
well as the structure and dimensionality of space-time HTJ. Microscopic black hole properties 
are reviewed in more detail in l|T5l [T6l. 

The search for black holes is based on i/s = 7 TeV pp collision data recorded by the Compact 
Muon Solenoid (CMS) detector at the LHC between March and October 2010, which corre- 
spond to an integrated luminosity of 34.7 ± 3.8 pb^^. A detailed description of the CMS ex- 
periment can be found elsewhere IITtI . The central feature of the CMS detector is the 3.8 T 
superconducting solenoid enclosing the silicon pixel and strip tracker, the electromagnetic cal- 
orimeter (ECAL), and the brass-scintilla tor hadronic calorimeter (HCAL). For triggering pur- 
poses and to facilitate jet reconstruction, the calorimeter cells are grouped in projective towers, 
of granularity Arj x Acp = 0.087 x 0.087 at central rapidities and 0.175 x 0.175 in the forward 
region. Here, the pseudorapidity rj is defined as — In (tan ^), where 6 is the polar angle with 
respect to the direction of the counterclockwise beam, and cp is the azimuthal angle. Muons are 
measured in the pseudorapidity window \r]\ < 2.4 in gaseous detectors embedded in the steel 
return yoke. 

The CMS trigger system consists of two levels. The first level (LI), composed of custom hard- 
ware, uses information from the calorimeters and muon detectors to select the most interesting 
events for more refined selection and analysis at a rate of up to 80 kHz. The software-based 
High Level Trigger (HLT) further decreases the rate to a maximum of ~ 300 Hz for data storage. 
The instantaneous luminosity is measured using information from forward hadronic calorime- 
ters IIT8I . 

We use data collected with a dedicated trigger on the total jet activity, Hj, where Hj is defined 
as the scalar sum of the transverse energies Ej of the jets above a preprogrammed threshold. 
At LI this jet threshold was 10 GeV, and the Hj threshold was 50 GeV. At HLT, the jet Ej 
threshold varied between 20 and 30 GeV, and the Hj threshold between 100 and 200 GeV. The 
trigger is fully efficient for the offline analysis selections described below. Energetic electrons 
and photons are also reconstructed as jets at the trigger level and are thus included in the Hj 
sum. 

Jets are reconstructed using energy deposits in the HCAL and ECAL, clustered using a collinear 
and infrared safe anti-A;^ algorithm with a distance parameter of 0.5 [19J. The jet energy reso- 
lution is AE/ E 100 %/ a/ E [GeV] © 5 %. Jets are required to pass quality requirements to re- 
move those consistent with calorimeter noise. Jet energies are corrected for the non-uniformity 
and non-linearity of the calorimeter response, as derived using Monte Carlo (MC) samples and 
collision data |20|. Jets are required to have Ej > 20 GeV before the jet-energy-scale correc- 
tions and to have < 2.6. Missing transverse energy K^is reconstructed as the negative of the 
vector sum of transverse energies in the individual calorimeter towers. This quantity is further 
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corrected to account for muons in the event, which deposit little energy in the calorimeters, and 
for the jet energy scale l|2T) . 

Electrons and photons are identified as isolated energy deposits in the ECAL, with a shape 
consistent with that expected for electromagnetic showers. Photons are required to have no 
matching hits in the inner pixel detector layers, while electrons are required to have a matching 
track. Electrons and photons are required to have Ej > 20 GeV and to be reconstructed in 
the fiducial volume of the barrel {\f]\ < 1.44) or the endcap (1.56 < \r]\ < 2.4). The ECAL 
has an ultimate energy resolution better than 0.5% for unconverted photons or electrons with 
transverse energies above 100 GeV [22 J. In 2010 collision data, for Ej > 20 GeV, this resolution 
is better than 1% in the barrel. 

Muons are required to have matched tracks in the central tracker and the muon spectrometer, 
to be within \tj\ < 2.1, be consistent with the interaction vertex to suppress backgrounds from 
cosmic ray muons, be isolated from other tracks, and have transverse momentum pj above 
20 GeV. The combined fit using tracks measured in the central tracker and the muon spectrom- 
eter results in pj resolution between 1% and 5% for pj values up to 1 TeV 

The separation between any two objects (jet, lepton, or photon) is required to be 



Black hole signal events are simulated using the parton-level BlackMax [231 generator (v2.01.03), 
followed by a parton-showering fragmentation with PYTHI A [24J (v6.420), and a fast parametric 
simulation of the CMS detector response [25J, which has been extensively validated for signal 
events using detailed detector simulation via GEANT4 l|26l . 

Several additional parameters govern black hole production and decay in the ADD model in 
addition to Md and n. For each value of Md, we consider a range of the minimum black hole 
masses, MJ^, between Md and the kinematic limit of the LHC. We assume that no parton- 
collision energy is lost in gravitational shock waves, i.e. it is all trapped within the event hori- 
zon of the forming black hole. We consider both rotating and non-rotating black holes in this 
analysis. Graviton radiation by the black hole is not considered. For most of the signal samples 
we assume full Hawking evaporation without a stable non-interacting remnant. 

The parameters used in the simulations are listed in Tablejljfor a number of characteristic model 
points. The MSTW20081o68 [27J parton distribution functions (PDF) were used. In addition we 
compare the BlackMax results with those of the CHARYBDIS 2 MC generator (vl.0.3) Il28ll29l . 
The two generators yield different values of total cross section, as BlackMax introduces an 
additional n-dependent factor applied on top of the geometrical cross section. The CHARYB- 
DIS cross sections are a factor of 1.36, 1.59, and 1.78 smaller than those from BlackMax for 
n = 2, 4, and 6, respectively. In addition, CHARYBDIS has been used to simulate black hole 
evaporation resulting in a stable non-interacting remnant with mass Md (this model is not im- 
plemented in BlackMax). In the generation, we use the Particle Data Group [5J definition of the 
Planck scale Mu- (Using another popular choice for Md from Dimopoulos and Landsberg [IT| 
would result in a suppression of the production cross section by a factor of 1.35, 5.21, or 9.29 
for n = 2, 4, or 6, respectively.) 

We employ a selection based on total transverse energy to separate black hole candidate events 
from the backgrounds. The variable Sj is defined as a scalar sum of the £7- of the N individual 
objects (jets, electrons, photons, and muons) passing the above selections. Only objects with 
Ej > 50 GeV are included in the calculation of Sj, in order to suppress the SM backgrounds 
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Table 1: Monte Carlo signal points for some of the model parameters probed, corresponding 
leading order cross sections (cr), and the minimum required values for the event multiplicity 

(N > N"^^) and Sj (S^"^), as well as the signal acceptance (A), the expected number of signal 
events (n^^^), the number of observed events (n'^^^'^) in data, the expected number of background 
events (n^^^), and the observed (cr^^) and expected (c^lxp.) limits on the signal cross section at 
95% confidence level. 
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and to be insensitive to jets from pile-up, while being fully efficient for black hole decays. 
Further, the missing transverse energy in the event is added to Sj, if the missing transverse 
energy value exceeds 50 GeV. Note that while £r is counted toward St, it is not considered in 
the determination of N. 

The main background to black hole signals arises from QCD multijet events. Other back- 
grounds from direct photon, W/Z+jets, and tt production were estimated from MC simu- 
lations, using the MadGraph |30| leading-order parton-level event generator with CTEQ6L 
PDF set [31 1, followed by pythia [24] parton showering and full CMS detector simulation via 
GEANT4 [26[. These additional backgrounds are negligible at large values of Sj and contribute 
less than 1% to the total background after the final selection. 

The dominant multijet background can only be estimated reliably from data. For QCD events, 
Sr is almost completely determined by the hard 2^2 parton scattering process. Further 
splitting of the jets due to final-state radiation, as well as jets due to initial-state radiation - most 
often nearly collinear with either incoming or outgoing partons - does not change the Sj value 
considerably. Consequently, the shape of the Sr distribution is expected to be independent of 
the event multiplicity N, as long as Sj is sufficiently above the turn-on region (i.e., much higher 
than N X 50 GeV). 

We confirmed the assumption of the St shape in variance of N using MC generators capable of 
simulating multijet final states from either matrix elements [32 [ or parton showers [24 J. This 
shape invariance offers a direct way of extracting the expected number of background events 
in the search for black hole production. The decay of black holes with a mass of a few TeV typi- 
cally results in events with multiplicity of half-a-dozen objects in the final state. The conjecture 
that the St shape is independent of the multiplicity has also been checked with data using the 
exclusive multiplicities of N = 2 and N = 3. Even in the presence of a signal, its contribution to 
these multiplicities is expected to be small and only seen at large values of St, so these samples 
can be used for the background prediction at higher multiplicities. Moreover, since dedicated 
analyses of the dijet invariant mass spectrum have been conducted ['T3'.'T4'[, we know that there 
are no appreciable contributions from new physics to the dijet final state up to invariant masses 
of about 1.5 TeV, which, for central jets, translates to a similar range of St- 

We fit the St distributions between 600 and 1100 GeV, where no black hole signal is expected, 
for data events with N = 2 and N = 3 using an ansatz function iol{x) > which is shown 
with the solid line in Fig. [l] To check the systematic uncertainty of the fit, we use two addi- 
tional ansatz functions, jz= — — ^rp- and ° flSl, which are shown as the upper and lower 

boundaries of the shaded band in Fig. [l] The default choice of the ansatz function was made 
based on the best-fit to the Sj distribution for N = 2. Additional systematic uncertainty arises 
from a slight difference between the best-fit shapes for N = 2 and N = 3. Nevertheless, the fits 
for these two exclusive multiplicities agree with each other within the uncertainties, demon- 
strating that the shape of the St distribution is independent of the final-state multiplicity. 

The St distributions for data events with multiplicities N > 3, 4, and 5 are shown in Figs.|2^, b, 
and c, respectively. The solid curves in the figures are the predicted background shapes, found 
by normalizing the fits of the N = 2 and 3 Sr distributions to the range St = 1000 — 1100 GeV, 
where no black hole signal contribution is expected. 

Since no excess is observed above the predicted background, we set limits on the black hole 
production. We assign a systematic uncertainty on the background estimate of 6% to 125% for 
the St range used in this search. This uncertainy comes from the normalization uncertainty (4 - 
12%, dominated by the statistics in the normalization region) added in quadrature to the uncer- 
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St (GeV) (GeV) 

Figure 1: Total transverse energy Sj, for events with the multiplicities of a) N = 2, and b) N = 3 
objects in the final state. Data are depicted as solid circles with error bars; the shaded band is 
the background prediction obtained from data (solid line) with its uncertainty. Non-multijet 
backgrounds are shown as colored histograms. Also shown is the predicted black hole signal 
for three different parameter sets. 

tainties arising from using various ansatz fit functions and the difference between the shapes 
obtained from the N = 2 and N = 3 samples. The integrated luminosity is measured with an 
uncertainty of 11% flS]. The uncertainty on the signal yield is dominated by the jet energy scale 
uncertainty of ~ 5% [20J which translates into a 5% uncertainty on the signal. An additional 
2% uncertainty on the signal acceptance comes from the variation of PDFs within the CTEQ6 
error set liMI . The particle identification efficiency does not affect the signal distribution, since 
an electron failing the identification requirements would be classified either as a photon or a 
jet; a photon failing the selection would become a jet; a rejected muon would contribute to the 
Bt- In any case the total value of Sj is not affected. 

We set limits on black hole production with the optimized Sj and N selections by counting 
events with Sj > S^™ and N > N'"'". We optimized the signal (S) significance in the presence 
of background (B) using the ratio S/ \/S + B for each set. The optimum choice of parameters is 
listed in Table [l| as well as the predicted number of background events, the expected number 
of signal events, and the observed number of events in data. Note that the background uncer- 
tainty, dominated by the choice of the fitting function, is highly correlated for various working 
points listed in Table[l]and also bin-to-bin for the Sj distributions shown in Figs.[l]and|2] 

We set upper limits on the black hole production cross section using the Bayesian method 
with flat signal prior and log-normal prior for integration over the nuisance parameters (back- 
ground, signal acceptance, luminosity) ||5l|33j. These upper limits at the 95% confidence level 
(CL) are shown in Fig. |3j as a function of M|^. For the three working points shown in the 
figure, the observed (expected) lower limits on the black hole mass are 3.5, 4.2 and 4.5 TeV (3.2, 
4.0, and 4.5 TeV), respectively. 

Translating these upper limits into lower limits on the parameters of the ADD model, we can 
exclude the production of black holes with minimum mass of 3.5 — 4.5 TeV for values of the 
multidimensional Planck scale up to 3.5 TeV at 95% CL. These limits, shown in Fig. |4| do not 
exhibit significant dependence on the details of the production and evaporation model. These 
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Figure 2: Total transverse energy Sj, for events with multiplicities a) N > 3, b) N > 4, and c) 
N > 5 objects in the final state. Data are depicted as solid circles with error bars; the shaded 
band is the background prediction (solid line) with its uncertainty. Also shown are black hole 
signals for three different parameter sets. 
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Figure 3: The 95% confidence level upper limits on the black hole production cross section 
(solid lines) and three theoretical predictions for the cross section (dotted lines), as a function 
of the black hole mass. 
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dimensional Planck scale Md for several benchmark scenarios. The area below each curve is 
excluded by this search. 



9 



are the first limits of a dedicated search for black hole production at hadron colliders. 
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Figure 5: Model-independent 95% confidence level upper limits on a signal cross section times 
acceptance for counting experiments with Sj > S^™ as a function of for (a) N > 3, (b) 
N > 4, and (c) N > 5. The blue (red) lines correspond to an observed (expected) limit for 
nominal signal acceptance uncertainty of 5%. 

Finally, we produce model-independent upper limits on the cross section times the acceptance 
for new physics production in high-Sr inclusive final states for N > 3, 4, and 5. Figure |5] shows 
95% CL upper limits from a counting experiment for Sj > S™" as a function of S j™, which can 
be used to test models of new physics that result in these final states. A few examples of such 
models are production of high-mass tt resonances |34| in the six-jet and lepton + jet final states, 
R-parity violating gluino decay into three jets, resulting in the six-jet final state 1.35 , and 
a class of models with strong dynamics, with a strongly produced resonance decaying into a 
pair of resonances further decaying into two jets each, resulting in the four-jet final state ||37 [ . In 
addition, these limits can be used to constrain black hole production for additional regions of 
the parameter space of the model, as well as set limits on the existence of string balls [38 j , which 
are quantum precursors of black holes predicted in certain string models. We have checked 
that for the black hole model parameters we probed with the dedicated optimized analysis, the 
sensitivity of the search in terms of the excluded black hole mass range exceeds that from the 
model-independent cross section limits by as little as 5 - 8%. Thus, model-independent limits 
can be used efficiently to constrain the allowed parameter space in an even broader variety of 
black hole models than we covered in this Letter. 

To conclude, we have performed the first dedicated search for microscopic black holes at a par- 
ticle accelerator and set limits on their production within a variety of models. The lower limits 
on the black hole mass at 95% CL range from 3.5 to 4.5 TeV for values of the Planck scale up 
to 3.5 TeV in the model with large extra dimensions in space. Additionally, we have produced 
model-independent limits on the production of energetic, high-multiplicity final states, which 
can be used to constrain a variety of models of new physics. 
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